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ABSTRACT 

We present in this paper the general real- and redshift-space clustering properties of galaxies as measured in the first data release 
of the VIPERS survey. VIPERS is a large redshift survey designed to probe the distant Universe and its large-scale structure at 
0.5 < z < 1.2. We describe in this analysis the global properties of the sample and discuss the survey completeness and associated 
corrections. This sample allows us to measure the galaxy clustering with an unprecedented accuracy at these redshifts. From the 
redshift-space distortions observed in the galaxy clustering pattern we provide a first measurement of the growth rate of structure at 
z - 0.8: /cTg = 0.47 ± 0.08. This is completely consistent with the predictions of standard cosmological models based on Einstein 
gravity. 

Key words. Cosmology; observations - Cosmology: large scale structure of Universe - Galaxies: high-redshift - Galaxies: statistics 
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i1. Introduction 

Galaxy redshift surveys have provided a wealth of information 
on the inhomogeneous universe, mapping the late -time develop- 
ment of the small metric fluctuations that existed at early times, 
and whose early properties can be viewed in the cosmic mi- 
crowave background (CMB). The growth of structure during this 
intervening period is sensitive both to the type and amount of 
dark matter, and also to the theory of gravity, so there is a strong 
motivation to make precise measurements of the rate of growth 
of cosmological structure (e.g. Jain & Khoury 2010). 

Of course, galaxy surveys do not image the mass fluctua- 
tions directly, unlike gravitational lensing. But the visible light 
distribution does have some advantages as a cosmological tool 
in comparison with lensing. The number density of galaxies is 
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vatory, Cerro Paranal, Chile, using the Very Large Telescope under 
programs 182.A-0886 and partly 070.A-9007. Also based on obser- 
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de la Recherche Scientifique (CNRS) of France, and the University of 
Hawaii. This work is based in part on data products produced at TER- 
APIX and the Canadian Astronomy Data Centre as part of the Canada- 
France-Hawaii Telescope Legacy Survey, a collaborative project of 
NRC and CNRS. The VIPERS web site is http://www.vipers.inaf.it/. 



sufliciently high that the density field of luminous matter can 
be measured with a finer spatial resolution, probing interesting 
nonhnear features of the clustering pattern with good signal-to- 
noise. The price to be paid for this is that the complicated biasing 
relation between visible and dark matter has to be confronted; 
but this is a positive factor in some ways, since understanding 
galaxy formation is one of the main questions in cosmology. 
Redshift surveys provide the key information needed to meet this 
challenge: global properties of the galaxy population and their 
variation with environment and with epoch. 

The final advantage of redshift surveys is that the radial in- 
formation depends on cosmological expansion and is corrupted 
by peculiar velocities. Although the lack of a simple method to 
recover true distances can be frustrating at times, it has come 
to be appreciated that this complication is in fact a good thing. 
The pecuUar velocities induce an anisotropy in the apparent clus- 
tering, from which the properties of the pecuUar velocities can 
be inferred much more precisely than in any attempt to measure 
them directly using distance estimators. The reason peculiar ve- 
locities are important is that they are related to the underlying 
linear fractional density perturbation 6 via the continuity equa- 
tion: 5 = -V • u, where u is the peculiar velocity field. This can 
be expressed more conveniently in terms of the dimensionless 
scale factor, ait), and the Hubble parameter, H(t), as 



V • u = -Hf6; f = 



dlnd 
d\na 



(1) 
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The growth rate can be approximated in most models by f(a) ^ 
^miay, where y ^ 0.545 in standard A-dominated models, but 
where models of non-standard gravity display a growth rate in 
which the effective value of y can differ by 30% (Linder & Cahn 
2007). 

The possibility of using the redshift-space distortion sig- 
nature as a probe of the growth rate of density fluctuations is 
one of the main reasons behind the recent burst of activity in 
galaxy redshift surveys. The first paper to emphasise this ap- 
plication as a test of gravity theories was the analysis of the 
VVDS survey by Guzzo et al. (2008), and subsequent work es- 
pecially by the SDSS LRG (Samushia et al. 2011), WiggleZ 
(Blake et al. 2011; Contreras et al. 2013), 6dFGS (Beutler et al. 
2012) and BOSS (Reid et al. 2012) surveys has exploited this 
method to make measurements of the growth rate at z < 1. 
The drive in this field is both towards higher accuracy (larger 
surveys) and a broader redshift range, and this provides the 
motivation for the VIMOS Public Extragalactic Redshift Sur- 
vey (http : //vipers . inaf . it). VIPERS is an on-going ESO 
Large Programme aimed at measuring redshifts for about 10^ 
galaxies at redshift 0.5 < z < 1.2. Its prime goal is an accu- 
rate measurement of the growth rate of large-scale structure at 
redshift one and beyond, but in general the survey is intended to 
provide robust and precise measurements of all properties of the 
galaxy population at an epoch when the Universe was about half 
its current age. 

This paper presents the initial analysis of the real-space 
galaxy clustering and redshift-space distortions in VIPERS, to- 
gether with the resulting impUcations for the growth rate. The 
data are described in Section 2; Section 3 describes the survey 
selection effects; Section 4 describes our methods for estimating 
clustering, which are tested on simulations in Section 5; Section 
6 presents the real-space clustering results; Section 7 gives the 
redshift-space distortions results, and Section 8 summarises our 
results and concludes. 

Throughout this analysis, if not specified otherwise, 
we assume a fiducial ACDM cosmological model with 
(Qm,£2;fc, vw) = (0.25,0,-1) and a Hubble constant of Hq = 
100/ikms-i Mpc-^ 



2. Data 

The VIPERS galaxy target sample is selected from the optical 
photometric catalogues of the Canada-France-Hawaii Telescope 
Legacy Survey Wide (CFHTLS-Wide, Goranova et al. 2009). 
VIPERS covers 24 deg^ on the sky, divided over two areas within 
the Wl and W4 CFHTLS fields. Galaxies are selected to a hmit 
of < 22.5, applying a simple and robust gri colour pre- 
selection to efficiently remove galaxies at z < 0.5. Coupled with 
a highly optimized observing strategy (Scodeggio et al. 2009), 
this allows us to double the galaxy sampUng rate in the redshift 
range of interest, with respect to a pure magnitude-Unoited sam- 
ple. At the same time, the area and depth of the survey result in 
a relatively large volume, 5 x 10^ h^^ Mpc^, analogous to that 
of the Two Degree Field Galaxy Redshift Survey (2dFGRS) at 
z ^ 0.1 (CoUess et al. 2001, 2003). Such a combination of sam- 
phng rate and depth is unique amongst current redshift surveys at 
z > 0.5. VIPERS spectra are collected with the VIMOS multi- 
object spectrograph (Le Fevre et al. 2003) at moderate resolu- 
tion (R = 210) using the LR Red grism, providing a wavelength 
coverage of 5500-9500A and a typical radial velocity error of 
cr, = 175(1 + z)kms-'. The full VIPERS area of 24 deg^ will 
be covered through a mosaic of 288 VIMOS pointings (192 in 



the Wl area, and 96 in the W4 area). A discussion of the sur- 
vey data reduction and management infrastructure is presented 
in Garilli et al. (2012). An early subset of the spectra used here 
is analysed and classified through a Principal Component Anal- 
ysis (PCA) in Marchetti et al. (2013). A complete description of 
the survey construction, from the definition of the target sample 
to the actual spectra and redshift measurements, is given in the 
parallel survey description paper (Guzzo et al. 2013). 

The data set used in this and the other papers of the early 
science release, will represent the VIPERS Pubhc Data Release 
1 (PDR-1) catalogue. It will be pubUcly available in the fall of 
2013. This catalogue includes 55358 redshifts (27935 in Wl and 
27423 in W4) and corresponds to the reduced data frozen in the 
VIPERS database at the end of the 2011/2012 observing cam- 
paign; this represents 64% of the final survey in terms of cov- 
ered area. A quaUty flag has been assigned to each object in the 
process of determining their redshift from the spectrum, which 
quantifies the reliability of the measured redshifts. In this analy- 
sis, we use only galaxies with flags 2 to 9 inclusive, correspond- 
ing to a sample with a redshift confirmation rate of 98%. The 
redshift confirmation rate and redshift accuracy have been esti- 
mated using repeated spectroscopic observations in the VIPERS 
fields (see Guzzo et al. 2013, for details). The catalogue, which 
we will refer to just as the VIPERS sample in the following, cor- 
responds to a sub-sample of 4587 1 galaxies with reUable redshift 
measurements. 

The redshift distribution of the sample is presented in Fig. 1 . 
We can see in this figure that the survey colour selection allows 
an efficient removal of galaxies below z = 0.5. It is important to 
notice that the colour selection does not introduce a sharp cut in 
redshift but a redshift window function which has a smooth tran- 
sition from zero to one in the redshift range 0.4 < z < 0.6, with 
respect to the full population of i' < 22.5 galaxies. This effect 
on the radial selection of the survey, which we refer to as the 
Colour Sampling Rate (CSR) in the following, is only present 
below z - 0.6. Above this redshift, the colour selection has no 
impact on the redshift selection and the sample becomes purely 
magnitude-limited at i' < 22.5 (Guzzo et al. 2013). If we weight 
the raw redshift distribution by the global survey completeness 
function described in the next sections, one obtains the N(z) rep- 
resented by the empty histogram in Fig. 1. For convenience, 
we scaled down the corrected N{z) by 40%, the average effective 
survey sampling rate, to aid the comparison between the shapes 
of the two distributions. The difference in shape between these 
two N(z) shows the effect of incompleteness in the survey, which 
is only significant at about z > 0.9 (see also Davidzon et al. 
2013). 

The observed redshift distribution in the sample can be well 
described by a function of the form 



iV(z)=A^--j exp[^-^^| \CSRiz), (2) 

in units of deg"^ ■ (Az = 0.03)"' and where (A,ZQ,a,(T} = 
(3.103, 0.191, 8.603, 1.448). The CSR is flie incompleteness in- 
troduced by the VIPERS colour selection. It is primarily a func- 
tion of redshift and can be estimated from the ratio between the 
number of galaxies with /' < 22.5 satisfying the VIPERS colour 
selection and the total number of galaxies with /' < 22.5 as a 
function of redshift. We cahbrated this function using the VLT- 
VIMOS Deep Survey Wide spectroscopic sample (VVDS- Wide, 
Garilli et al. 2008) which has a CFHTLS-based photometric cov- 
erage and depth that is similar to that of VIPERS, but which is 
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Fig. 1. Redshift distribution of the combined W1+W4 galaxy sample 
when including only reliable redshifts (filled histogram) and that cor- 
rected for the full survey completeness (empty histogram). The curve 
shows the best-fitting template redshift distribution given by Eq. 2 ap- 
plied to the uncorrected observed distribution. 



free from any colour selection (see Guzzo et al. 2013, for de- 
tails). 

The fitting of N{z) is important in measuring galaxy cluster- 
ing: the form of the mean redshift distribution must be followed 
accurately, but features from large-scale structure must not be al- 
lowed to bias the result. We discuss this issue in detail in Section 
5. 



3. Angular completeness 

3.1. Slit assignment and footprint 

To obtain a sample of several square degrees with VIMOS, one 
needs to perform a series of individual observations or point- 
ings. The VIPERS strategy consists in covering the survey area 
with only one pass. This has been done in order to maximise 
the volume probed. The survey strategy and the fact that the VI- 
MOS field-of-view is composed of four quadrants delimited by 
an empty cross, create a particular footprint on the sky which 
is reproduced in Figs 4 and 5. In each pointing, slits are as- 
signed to a number of potential targets which meet the survey 
selection criteria. This is shown in Fig. 2, which illustrates how 
the slits are positioned in the pointing W1P082. Given the sur- 
face density of the targeted population, the multiplex capability 
of VIMOS, and the survey strategy, a fraction of about 45% of 
the parent photometric sample can be assigned to slits. We de- 
fine the fraction of target which have a measured spectrum as 
the Target Sampling Rate {TSR) and the fraction of observed 
spectra with reliable redshift measurement as the Spectroscopic 
Sampling Rate (SSR). The number of slits assigned per pointing 
is maximized by the SSPOC algorithm (Bottini et al. 2005), but 
the elongated size of the spectra means that the resulting sam- 
pling rate is not uniform inside the quadrants. The dispersion 
direction of the spectra in VIPERS are aligned with the Dec di- 
rection and consequently, the density of spectra along this direc- 
tion is lower with respect to that along the RA direction. This 



particular sampling introduces an observed anisotropic distribu- 
tion of pair separation, which has to be accounted for in order to 
measure galaxy clustering correctly. 

The two empty stripes between the four quadrants in each 
pointing introduce a particular pattern in the measured correla- 
tion functions if not accounted for. We correct for that by ap- 
plying detailed binary masks of the spectroscopic observations 
to a random sample of unclustered objects, so that both data and 
random catalogues contain no objects in these stripes. These 
masks account for the detailed VIMOS field-of-view geometry 
as well as for the presence of vignetted areas at the boundaries 
of the pointings. On top of these spectroscopic masks, we apply 
a set of photometric masks which discard areas where the par- 
ent photometry is affected by defects such as large stellar haloes 
and where the survey selection is compromised (see Guzzo et al. 
2013). 

3.2. Small-scale incompleteness 

We can characterise the amount of missing small-scale angular 
pairs induced by the VIPERS spectroscopic strategy, by measur- 
ing the angular pair completeness as a function of angular sep- 
aration. This quantity, defined as the ratio between the number 
of pairs in the spectroscopic sample and that in the parent pho- 
tometric sample, can be written in terms of angular two-point 
coiTelation functions as (Hawkins et al. 2003) 



(3) 



1 _ 1 -H W,(0) 

wHff) ~ 1 + Wpi9) ' 

where Ws{0) and Wp(0) are respectively the angular correlation 
function of the spectroscopic and parent samples. This function 
is shown in Fig. 3. No significant difference is seen between the 
Wl and W4 fields, as expected. The amount of missing angular 
pairs is only significant below = 0.03 deg, which corresponds 
to a transverse comoving scale of about 1 /i"' Mpc at z = 0.8. 

This fraction varies with redshift, although in practice we 
cannot measure it at different redshifts since we do not have a 
measured redshift for all galaxies in the parent sample. For this 
reason we use the global w^(ff) (averaged over all observed red- 
shifts) to correct for the small-scale angular incompleteness ef- 
fect. We will show in Section 5 that the level of systematic error 
introduced by using w^{0) instead of w^{6\z) is very small, of 
the order of a few percent. When measuring the angular cor- 
relation functions, we include the completeness weights intro- 
duced in the following section, in a similar way as for the thee- 
dimensional coiTelation function estimation. 

It is important to mention that the small-scale angular incom- 
pleteness effect is a general issue for large galaxy redshift sur- 
veys, in which one has to deal with the mechanical constraints 
of multi-object spectrographs and survey strategy. The incom- 
pleteness due to slit assignment in VIPERS is to some extent 
similar to the fibre collision problem in surveys using fibre spec- 
troscopy such as 2dFGRS or SDSS, while the magnitude of the 
effect is much more severe in our case. Recently, a new method 
has been developed to accurately correct for fibre collision (Guo 
et al. 2012). Although this method is quite general, it is not appli- 
cable here. The exclusion between spectroscopically observed 
objects in VIPERS is essentially uni-directional, meaning that 
not all close pairs are excluded. Therefore calculations such as 
that shown in Fig. 4 are possible from the set of one-pass ob- 
servations, whereas they can only be included for SDSS where 
overlapping observations are included. Thus we need to revise 
the coiTection methods developed for such surveys to apply them 
to VIPERS. 
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Fig. 2. Illustration of the slit assignment in pointing W1P082. The slits are shown in red and associated rectangles represent the typical dispersion 
of the spectra. All objects meeting the survey selection criteria (potential spectroscopic targets) are represented by black circles. 



3.3. Large-scale incompleteneness 

In addition to the non-uniform sampling inside the pointings, the 
survey has variations of completeness from quadrant to quad- 
rant. This incompleteness is the combined effect of the TSR and 
SSR. The latter, which characterises our ability of determining 
a redshift from a galaxy spectrum, is determined empirically as 
the ratio between the number of reliable redshifts and the total 
number of observed spectra. The TSR and SSR in each quad- 
rant is shown in in Figs 4 and 5. From these figures one can 
see clearly that both TSR and SSR functions vary according to 
the position on the sky, although SSR tends to have stronger 
variations. The variations of TSR reflect the changes in angu- 
lar galaxy density in the parent catalogue. Indeed, because of 
the finite maximum number of slits that can be assigned and the 
fact that each quadrant has a different number of potential tar- 
get, the less dense quadrants tend to be better sampled than the 
denser ones. On the other hand, variations in observational con- 
ditions from pointing to pointing induce changes in SSR. These 
different observational conditions translate into variations of the 
signal-to-noise of the measured spectra and so in our ability of 
extracting a redshift measurement from them. These effects are 
taken into account in the clustering estimation by weighting each 
galaxy according to the reciprocal of the TSR and SSR. 



4. Clustering estimation 

We characterise the galaxy clustering in the VIPERS sample by 
measuring the two-point statistics of the spatial distribution of 
galaxies in configuration space. We estimate the two-point cor- 
relation function ^(r) using the Landy & Szalay (1993) estimator 



GG{r) - IGRir) + RR{r) 
RR(r) ' 



(4) 



where GG(r), GR(r), and RR(r) are respectively the normalized 
galaxy-galaxy, galaxy-random, and random-random number of 
pairs with separation inside [r - Ar/2, r + Ar/2]. Note that here 
r is a general three-dimensional galaxy separation, not specif- 
ically the real-space separation. This estimator minimises the 
estimation variance and circumvent discreteness and finite vol- 
ume effects (Landy & Szalay 1993; Hamilton 1993). A random 
catalogue must be constructed in this estimator, whose aim is to 
accurately estimate the number density of objects in the sample. 
It must be an unclustered population of objects with the same 
radial and angular selection functions as the data. In this analy- 
sis, we use random samples with 20 times more objects than in 
the data to minimise the shot noise contribution in the estimated 
correlation functions. 
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Fig. 4. Variations of the target success rate (TSR) with quadrants. The TSR quantifies our ability of obtaining spectra from the potential 
targets meeting the survey selection in the parent photometric sample. The quadrants filled in black correspond to failed observations where no 
spectroscopy has been taken. 
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Fig. 5. Variations of the spectroscopic success rate (SSR) with quadrants. The SSR quantifies our ability of determining galaxy redshifts from 
observed spectra. The quadrants filled in black correspond to failed observations where no spectroscopy has been taken. 



VIPERS has a complex angular selection function which has 
to be taken into account carefully when estimating the correla- 
tion function. For this, we weight each galaxy by the survey 
completeness weight, as well as each pair by the angular pair 
weights described in the previous section (Eq. 3). The survey 
completeness weights correspond to the inverse of the effective 
sampling rate ESR in each quadrant Q, defined as 

w(Q) = ESR HQ) = iSSRiQ) X TSRiQ)) '. (5) 

By applying these weights we effectively up-weight galaxies in 
the pair counts. It is important to note that here we keep the spa- 
tial distribution of the random objects uniform across the survey 
volume. We recall that survey completeness weights account for 
the quadrant-to-quadrant variations of the survey completeness 
described in Section 3.3 but do not correct for the internal quad- 
rant incompleteness. For that we use the angular pair weights 
vA(ff) which are applied to the GG pair counts. In principle the 
ESR is also a function of redshift and galaxy type (see Davidzon 
et al. 2013). But given the statistics of the sample it is impossible 
to measure the additional dependence of this function on redshift 
and galaxy properties. Therefore, we decided to only account for 



its quadrant-to-quadrant variations. We discuss the accuracy of 
this approximation in Section 5. 

Additional biases can arise if the radial selection function 
exhibits strong variations with redshift. The effect is partic- 
ularly significant for magnitude-limited catalogues covering a 
large range of redshifts and in which the radial selection func- 
tion rapidly drops at high redshift. In that case, the pair counts 
is dominated by nearby, more numerous objects: distant objects, 
although probing larger volumes, will have less weight. To ac- 
count for this we use the minimum variance estimator of Davis 
& Huchra (1982) for which the galaxy counts are essentially 
weighted by the inverse of the volume probed by each galaxy. 
This weighting scheme, usually referred as the Jt, weighting, is 
defined as (Hamilton 1993) 

w'\z,s)^ ..,1, . , , , (6) 
1 + n{z)AnJi{s) 

where z is the redshift of the object, s is the redshift-space pair 
separation, h{r) the galaxy number density at z and Jt,{s) is de- 
fined as 

Us)^ r s'^^{s')ds' . (7) 
Jo 
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Fig. 3. Completeness fraction of angular galaxy pairs due to the 
slit- spectroscopy strategy in the Wl and W4 fields for all galaxies at 
0.5 < z < 1.0. This has been obtained from the parent and spectroscopic 
sample angular correlation function. 



Each pair is then weighted by, 



(8) 



However, we find that applying 73 weighting does not signifi- 
cantly change the ampHtude and shape of the correlation func- 
tion in our sample, and tends to produce noisy correlation func- 
tions especially for high-redshift sub-samples. We thus decided 
not to apply this correction in this analysis. 

The final weight assigned to GG, GR, and RR pairs combine 
the survey completeness and angular pair weights as 



No Na 



G/^w^J^wKe,)©,,^ 

'=1 j=l 

A-,, Ng 



(9) 
(10) 
(11) 



where 0,j(r) is equal to unity for r,j in [r - Ar/2, r + Ar/2] and 
null otherwise. 

We measure correlation functions using both linear and log- 
arithmic binning. We define the separation associated with each 
bin as the bin centre and as the mean pair separation inside the 
bin, respectively for the linear and logarithmic binning (Zehavi 
et al. 201 1). The latter definition is more accurate than using the 
bin centre, in particular at large r when the bin size is large. 

The galaxy real-space correlation function ^(r) is not directly 
measurable from redshift survey catalogues because of galaxy 
peculiar velocities that affect redshift measurements. Galaxy 
peculiar velocities introduce distortions in the galaxy clustering 
pattern and as a consequence we can only measure redshift-space 



quantities. We measure the anisotropic redshift-space correla- 
tion function ^{rp,ji) in which the redshift-space galaxy separa- 
tion vector has been divided in two components, and tt, re- 
spectively perpendicular and parallel to the line-of-sight (Fisher 
et al. 1994). This decomposition, which assumes the plane- 
parallel approximation, allows us to isolate the effect of pecu- 
liar velocities as these modify only the component parallel to 
the line-of-sight. The redshift-space distortions can then be mit- 
igated by integrating ^{rp,n) over the tt, thus defining the pro- 
jected correlation function 



-r 

J -Tim 



Wp(rp) using an optimal value of 7r„ 



(12) 



We measure 

40/!"' Mpc, allowing us to reduce the underestimation of the 
amplitude of Wp(rp) on large scales without including too much 
noise. The projected correlation function allows us to measure 
real-space clustering (but see the later parts of Section 5). Fi- 
nally, to combine the correlation function measurements coming 
from the two fields, we measure the number-density weighted 
mean of one plus the correlation functions ^{rp,n) in Wl and 
W4. 



5. Tests of the clustering estimation 

5.1. Simulation data 

To test the robustness of our clustering estimation we make use 
of a large number of mock galaxy samples, which are designed 
to be a realistic match to the VIPERS sample. We create two 
sets of mock samples based on the Halo Occupation Distribu- 
tion (HOD) technique. These two sets only differ by the input 
halo catalogue that has been used. In the first set of mocks, we 
used the haloes extracted from the MultiDark dark matter N- 
body simulation (Prada et al. 2012). This simulation, which as- 
sumes a flat ACDM cosmology with (f2,„, Qa, O^, h, n, erg) = 
(0.27, 0.73, 0.0469, 0.7, 0.95, 0.82), covers a volume of 
1 h^^ Gpc^ using = 2048"^ particles. In the simulation, the 
haloes have been identified using a friends -of-friends algorithm 
with a relative linking length of b = 0.17 times the interparti- 
cle separation (i.e. 0.083/;"' Mpc) . The mass limit to which 
halo catalogues are complete is lO"'^ /i"' Mg. Because this lim- 
iting mass is too large to host the faintest galaxies observed with 
VIPERS, we use the method of de la Torre & Peacock (2012) 
to reconstruct haloes below the resolution limit. This method 
is based on stochastically resampling the halo number density 
field using constraints from the conditional halo mass function. 
For this, one needs to assume the shapes of the halo bias fac- 
tor and halo mass function at masses below the resolution limit 
and use the analytical formulae obtained by Tinker et al. (2008, 
2010). With this method we are able to populate the simulation 
with low-mass haloes with a sufficient accuracy to have unbi- 
ased galaxy two-point statistics in the simulated catalogues (see 
de la Torre & Peacock 2012, for details). The minimum re- 
constructed halo mass we consider for the purpose of creating 
VIPERS mocks is 10'" /i"' Mq. 

We then apply to the complete halo catalogues the algorithm 
presented in Carlson & White (2010) to remap halo positions 
and velocities in the initial simulation cube onto a cuboid of the 
same volume but different geometry. This is done to accom- 
modate a maximum number of disjoint VIPERS Wl and W4 
fields within the 1 h^^ Gpc^ volume of the simulation. This pro- 
cess allows us to create 26 and 31 independent lightcones for 
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Wl and W4 respectively over the redshift range 0.4 < z < 1.3. 
The lightcones are built by considering haloes from the different 
snapshots, disposing them according to their distance from the 
coordinate origin of the lightcone. The lightcones are then pop- 
ulated with galaxies using the HOD technique. In this process, 
we populate each halo with galaxies according to its mass, the 
mean number of galaxies in a halo of a given mass being given 
by the HOD. It is common usage to differentiate between central 
and satellite galaxies in haloes. While the former are put at rest at 
halo centres, the latter are randomly distributed within each halo 
according to a NFW radial profile. The halo occupation func- 
tion and its dependence on redshift and luminosity/stellar mass 
must be precisely chosen in order to obtain mock catalogues with 
realistic galaxy clustering properties. We calibrated the halo oc- 
cupation function directly on the VIPERS data. We performed 
an analytic HOD modelling of the projected correlation func- 
tion for different samples selected in luminosity and redshift that 
we will present in Section 6. We obtain from this a series of 
HOD parameters at different redshifts and for different cuts in B- 
band absolute magnitude, which we then interpolate to obtain a 
general redshift- and B-band absolute magnitude-dependent halo 
occupation function {Ng.^j(m\z, Mb))- We use the latter function 
to populate the haloes with galaxies. Finally, we add veloci- 
ties to the galaxies and measure their redshift-space positions. 
While the central galaxies are assigned the velocity of their host 
halo, satellite galaxies have an additional random component for 
which each Cartesian velocity component is drawn from a Gaus- 
sian distribution with a standard deviation that depends on the 
mass of the host halo. Details about the galaxy mock catalogue 
construction are given in Appendix A. 

The second set of mocks that we built is based on halo cat- 
alogues created with the Pinocchio code' (Monaco et al. 2002). 
This code follows the evolution of a set of particles on a reg- 
ular grid using an elUpsoidal model to compute collapse times 
and identify dark matter halos, and the Zel'dovich approxima- 
tion to displace the haloes from their initial position. While the 
recovery of haloes works well on an object-by-object basis, their 
positions and velocities on scales below 10 Mpc suffer by 
the lack of accuracy of the Zel'dovich approximation. The halo 
positions and velocities obtained with this method are less accu- 
rate than those from the N-body simulation, and the halo cluster- 
ing is generally underestimated on scales below 3 Mpc (e.g. 
Monaco et al. 2002). However this approach has the advantage 
of being very fast and can be used to generate a large number of 
independent halo catalogue realizations. We created 200 inde- 
pendent halo mock realizations assuming the same cosmology 
as the MultiDark N-body simulation. The remaining steps in 
generating galaxy mock samples are similar to those used for the 
mocks based on the MultiDark simulation. The only difference 
is that here we do not need to divide each simulation into sub- 
volumes to generate different lightcones: we can directly create 
volumes of the size of the lightcones. 

The final step in obtaining fully realistic VIPERS mocks is 
to add the detailed survey selection function. The procedure 
that we follow is similar to that used in the VVDS and zCOS- 
MOS surveys, which were also based on VIMOS observations 
(Meneux et al. 2006; lovino et al. 2010; de la Torre et al. 2011). 
We start by applying the magnitude cut /' < 22.5 and the effect 
of the colour selection on the radial distribution of the mocks. 
The latter is done by depleting the mocks at z < 0.6 so as to re- 



We have used in this analyis a new version of this code, optimised 
to work on massively parallel computers, which is described in Monaco 
etal. (2013). 



produce the CSR. The mock catalogues that we obtain are then 
similar to the parent photometric sample in the data. We next 
apply the slit-positioning algorithm with the same setting as for 
the data. This allows us to reproduce the VIPERS footprint on 
the sky, the small-scale angular incompleteness and the varia- 
tion of TSR across the fields. Finally, we deplete each quadrant 
to reproduce the effect of the SSR. Thus we are able to produce 
realistic mock galaxy catalogues that contain the detailed survey 
completeness function and observational biases of VIPERS. 

5.2. Effects of systematics on the correlation function 

5.2.1 . Effects related to the radial selection function 

We first study the impact on our correlation function measure- 
ments of using different methods to estimate the radial selec- 
tion. A key aspect in three-dimensional clustering estimation is 
to have a smooth and unbiased redshift distribution from which 
the random sample can be drawn. In particular, when the data 
sample used to estimate the radial distribution is not very large, 
one generally has to deal with strong features associated with 
prominent structures; these must not be allowed to induce spuri- 
ous clustering in the random sample. 

There are several empirical methods for avoiding this prob- 
lem. One can for instance interpolate the biimed observed distri- 
bution using cubic splines, filter the observed distribution with 
a kernel sufficiently large to erase the strong features in the dis- 
tribution, or fit the observed distribution with a smooth template 
N{z) and then randomly sample it. In general most of the meth- 
ods are parametric and have to be calibrated. An alternative non- 
parametric method is the Vmax method. This method consists in 
randomly sampling the maximum volumes V^ax probed by each 
galaxy in the survey (e.g. Kovac et al. 2010; Cole 2011). The 
Vmax value for each galaxy corresponds to the volume between 
the minimum and the maximum redshifts Zmin and Zmax at which 
the galaxy is observable in the survey. 

Fig. 6 applies three such approaches to estimate the galaxy 
radial distribution in the combined W1-I-W4 sample: the analyt- 
ical N(z) of Eq. 2; the Gaussian filtering method; and the Vmax 
method. This figure shows the recovered N(Dc) in the random 
sample with each method, with being the radial comoving 
distance; in practice we work with N(Dc) instead of N(z). We 
find that the methods give different estimates of the radial distri- 
bution. In the case of the the Gaussian filtering, a kernel size of 
150/;"' Mpc is needed to smear out the peaks in the distribution, 
otherwise the recovered NiDc) is still affected by large struc- 
tures in the field - particularly by that at D^. ^ 1600 Mpc. As 
expected, the filtering method tends to artificially broaden the 
N(Dc) distribution, whereas the analytical and Vmax methods are 
much smoother by construction and do not broaden the N(Dc). 
We find that the Vmax estimate shows a shghtly flatter distribu- 
tion at the level of the peak of the distribution, which seems vi- 
sually to be more consistent with the data. In Fig. 7 we show 
the effect of using these different estimates of the radial distribu- 
tion on the shape of the measured correlation function. Gaussian 
filtering with a kemel size of 150 /i ' Mpc and analytical Niz) 
estimates both yield sUghtly smaller ampUtudes of the projected 
correlation function on scales of above 10 h^^ Mpc than the Vmax 
method. Gaussian filtering with a kernel size of 100 /z"' Mpc 
globally underestimates the clustering amplitude on Wp{rp) as 
expected, by about 5%. The analytical and Vmax methods give 
very similar answers, except on scales above 5 Mpc where 
the former tends to produce a smaller clustering amplitude by 
5 - 15% with respect to the latter. This comparison shows that 
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Fig. 6. Comparison of different estimators of the radial distribution in 
the combined W1+W4 sample. The filled histogram shows the number 
of galaxies in fine bins of radial comoving distance. The different curves 
correspond to random radial distribution realizations normalized to the 
number of objects in the data, obtained using the Vn,;,x (solid), Gaussian 
filtering (dashed and dot-dashed), or analytical (dotted) methods. The 
vertical line shows the minimum redshift considered in this analysis, i.e. 
z = 0.5. 



the Vmax method is more robust as it uniquely allows us to re- 
store some correlation signal at large separation. For this reason 
and the fact that it is non-parametric we finally decided to use 
the Vtnax estimate to measure two-point correlation functions. 



5.2.2. Effects related to the angular selection function 

The most crucial aspect of the galaxy clustering estimation in 
VIPERS is to account for the angular selection function. We 
test our methodology and the different assumptions discussed 
in Section 4 using the MultiDark mock samples. We measure 
the accuracy with which we can estimate the two-point correla- 
tion function, by confronting the two-point correlation functions 
measured in the parent catalogues with those measured in the 
observed mocks when different completeness corrections are in- 
cluded. We measure the average relative difference between the 
corrected observed mocks and the parent measurement for dif- 
ferent statistics. For this test, we consider two galaxy samples 
encompassing respectively all galaxies in the redshift intervals 
0.5 < z < 0.75 and 0.75 < z < 1.0, using the same redshift dis- 
tribution in the parent and observed mock samples to construct 
the radial selection function of the random sample. 

It is common usage in clustering analysis to account for 
the angular survey completeness by down-weighting the ran- 
dom pair counts. This is usually done by keeping the galaxy 
counts unweighted and depleting the random sample so as to re- 
produce the survey angular completeness. The same effect can 
be achieved by using a uniform angular distribution of random 
objects but weighting each of them by the inverse of the weight 
defined in Eq. 5. If we do that and set all the angular pair weights 
to unity, we obtain the systematic error on Wp(rp) shown with the 
dotted curves in Fig. 8. We concentrate first on the results in the 
interval 0.5 < z < 0.75. We can see in this figure that the recov- 
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Fig. 7. Impact of the use of different estimators of the radial distri- 
bution on the shape of the projected correlation function. The projected 
correlation functions obtained using the V„,a;t (solid), Gaussian filtering 
(dashed and dot-dashed), or analytical (dotted) methods are shown in 
the top panel, while the relative fractional differences with respect to 
the Kmax method are presented in the bottom panel. 



ered clustering with this method is underestimated by about 10% 
at about I < rp < 20 Mpc, and then drops rapidly to 35% be- 
low. The strong underestimation on small scales is due to the 
small-scale angular incompleteness effect inside the quadrants. 
The approach of modulating the random density is dubious in 
the context of VIPERS, since it treats the sampling variations as 
a pattern imposed on the large-scale structure. But because of the 
VIMOS slit allocation, these variations are strongly coupled with 
the true clustering (i.e. the observed sky distribution of VIPERS 
galaxies is rather uniform). It is therefore safer if we keep the 
random sample uniform but upweight the galaxies as described 
in section 4. In this case, we obtain the dot-dashed lines in Fig. 
8: these represent an improved estimation of Wp{rp), reducing 
the underestimation by 5 - 6%. As expected, further including 
the angular pairs weights permits us to remedy in part the un- 
derestimation on scales below 1 /i"' Mpc, where the systematic 
error reaches 15% (solid lines). 

So far, we have used the global survey completeness and an- 
gular weights, i.e. neglecting the redshift dependence. As an ex- 
ercise we use the redshift information from the parent mocks to 
compute the true redshift-dependent weights and we obtain the 
dashed lines in the figure. Including the redshift dependence in 
the weights has the effect of improving the recovery of the pro- 
jected correlation function by about 2% over all probed scales. 
However, this improvement is rather modest - indicating that the 
use of the redshift-independent weights is a good approximation. 
Our best estimate of Wp{rp) therefore allows us to recover the 
true correlation function of the mocks at 0.5 < z < 0.75 with 
about 7% and 16% underestimation respectively above and be- 
low 1 /i ' Mpc. In the redshift interval 0.75 < z < 1 (shown in 
th bottom panel of Fig. 8), we find the same behaviour except 
that the correlation function is globally better recovered with an 
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Fig. 8. Systematic error on the projected correlation function and im- 
pact of different corrections. This is calculated considering all VIPERS 
galaxies in the redshift intervals 0.5 < z < 0.75 (top panel) and 
0.75 < z < 1 (bottom panel). 



underestimation smaller than 2 - 3% at > 0.6 h ' Mpc with 
the best method. 



This test demonstrates that our methodology gives an accu- 
rate estimate of the galaxy clustering in VIPERS, even if there 
remains some residual systematic errors of up to 7% on the 
scales above 1 Mpc and 15% on smaller scales. We find that 
the effect varies with redshift, being more important at the lowest 
redshifts probed by VIPERS. Overall these systematics remain 
within the sample variance errors, shown with shaded regions in 
Fig. 8. 



5.3. Impact of possible residual zero-point uncertainties in 
the photometry 

At the time of writing, photometry from the latest CFHTLS 
release (T0007) has become available (Hudelot et al. 2012). 
We have compared magnitudes and colours of objects in the 
VIPERS sample with the new CFHTLS-T0007 photometry. For 
VIPERS, the most important feature of T0007 compared to pre- 
vious releases is that each tile in the CFHTLS has now been 
rescaled to an absolute calibration provided by a new photomet- 
ric pre-survey taken at CFHT for this purpose. In addition, in or- 
der to ensure that seeing variations between tiles and filters are 
correctly accounted for, this scaling has been done using aper- 
ture fluxes that are rescaled based on the seeing on each indi- 
vidual tile; detailed tests at Terapix have shown that mag_auto 
magnitudes, which are affected by seeing variations, are not suf- 
ficiently precise for the percent-level photometric accuracy that 
is the objective of T0007. 

An important consequence of this work for VIPERS is that 
the effect of seeing variation and photometric calibration er- 
rors are now cleanly separated; the stellar-locus fitting technique 
used to define the VIPERS selection using colours based on 
mag_auto magnitudes mixes both these effects. To estimate the 
size of colour and magnitude offsets between T0007 and the ac- 
tual VIPERS selection (based on T0005) colours of stars on each 
VIPERS tiles measured from Terapix IQ20 magnitudes (used to 
calibrate T0007) and from mag_auto magnitudes in both releases 
have been compared. We find that these offsets shift the colour- 
colour locus we devised to remove lower-redshift z < 0.5 galax- 
ies (Guzzo et al. 2013). 

We test the effect of these possible variations of the colour 
selection across the fields in the context of galaxy clustering es- 
timation. For this we use photometric redshifts and quantify the 
variations in N{z) due to tile-to-tile variations of the colour se- 
lection, assuming the T0007 photometry as the reference. When 
comparing the N{Zphor) in the different tiles, we find that the red- 
shift distribution varies in shape and amplitude at z < 0.6 but 
only in amplitude above. The typical amplitude variations are 
of the order of about 5% (Guzzo et al. 2013). We then measure 
the ratio between the N(z) per tile and that averaged over the 
fields and use it as a redshift-dependent correction factor To test 
how these variations of the colour selection affect the measured 
correlation function, we vary the N{z) in the random sample for 
each quadrant using the correction factor previously defined on 
the averaged N(z). The projected correlations obtained with and 
without this correction are shown in Fig. 9. 

We can see that the correction has the effect of decreas- 
ing the amplitude of the projected correlation function by about 
2 - 4% on scales below 10 Mpc. We find a similar effect 
on the redshift-space angle-averaged correlation function ^{s). 
The amplitude and direction of the systematic effect follows our 
expectations, since spurious tile-to-tile fluctuations, if not prop- 
erly corrected, enhance the amplitude of clustering. This test 
suggests that indeed such tile-to-tile variations of colour selec- 
tion are present in the data. It is interesting to note that this 
systematic effect goes in the opposite direction to the effects of 
slit-positioning and associated incompleteness. In the end, be- 
cause this possible effect remains very small, we do not attempt 
to correct it for the clustering analysis. 



6. Real-space clustering 

Before studying redshift-space distortions in VIPERS, we begin 
by looking at the clustering in real space. The projected corre- 
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colour selection in the estimation of the projected correlation function, 
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lation function for all galaxies in the redshift range 0.5 < z < 1 
is shown in Fig. 10. It is measured in logarithmic bins of 
Alogrp - 0.2 over the scales 0.1 < rp < 30 h'^ Mpc. The 
error bars are estimated from the MultiDark mocks. 

The measured Wp(rp) functions in the Wl and W4 fields are 
very similar, in particular on scales below 5 /i"' Mpc. The com- 
bined projected correlation function in this redshift interval gives 
an accurate probe of the clustering up to scales of about 30/;"' 
Mpc. We can compare the galaxy projected correlation function 
to predictions for the mass non-linear correlation function and 
thus estimate the global effective linear bias of these galaxies. 
We use the HALOFIT (Smith et al. 2003) prescription for the 
non-linear mass power spectrum to compute the projected cor- 
relation function of mass at the mean redshift of the sample. By 
comparing the amplitudes of the measured galaxy and predicted 
mass correlations on scales of rp > 2/i"' Mpc, and assuming a 
linear biasing relation of the form obtain a 

linear bias of bi - 1.35 + 0.01. 

In order to make a detailed interpretation of the observed 
clustering of galaxies and produce realistic mock samples of the 
survey, we model our Wp(rp) measurements within the context 
of the HOD (Seljak 2000; Peacock & Smith 2000; Berlind & 
Weinberg 2002; Cooray & Sheth 2002). This method defines 
the mean distribution of galaxies within haloes; under the as- 
sumption of the abundance, large-scale bias, and density pro- 
file of haloes, one can then completely specify the clustering of 
galaxies and predict Wp(rp). We define four B-band absolute 
magnitude-threshold samples in the redshift bin 0.7 < z < 0.9 in 
which we measured Wp{rp). We model the projected correlation 
functions using HOD formalism, within a flat ACDM cosmol- 
ogy with parameters given by WMAP5 (Komatsu et al. 2009). 
We restrict the fit to scales above rp = 0.3/;"' Mpc and below 
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Fig. 10. Projected correlation functions of VIPERS galaxies in the 
redshift interval 0.5 < z < 1 for the individual Wl and W4 fields as 
well as for the combined sample. As a comparison, the ±lcr dispersion 
among the mean Wp(rp) in the moclcs is shown with the shaded region 
and the non-linear mass prediction in the assumed cosmology with the 
dotted curve. 
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Fig. 11. Measured and best-fitting HOD model projected correla- 
tion functions Wp(rp), for different luminosity-threshold subsamples of 
galaxies at 0.7 < z < 0.9. The Wp(rp) for the Mb - 5 log(/i) < -20.5, 
Mg-51og(/i) < -21.0, andMB-51og(/!) < -21.5 cases have been mul- 
tiplied by respectively 2, 3, and 4 to improve the clarity of the figure. 



rp - 30/;"' Mpc and correct empirically the measured projected 
correlation function for the residual underestimation at different 
scales, using the ratio between the parent and recovered Wp{rp) 
in the observed mocks for the same galaxy selection. We assume 
that there is negligible error in taking this small correction to be 
independent of cosmology. In the fitting procedure we used both 
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Fig. 12. The dependence of M^i„ and Mi HOD parameters on redshift 
and absolute magnitude threshold in S-band. The absolute magnitude 
threshold samples are plot in terms of their implied number density iig. 
The VIPERS results are compared to the pervious measurements per- 
formed in the DEEP2 (Zheng et al. 2007), BOSS (White et al. 2011), 
and CFHTLS (Coupon et al. 2012) surveys. 



the sample number density and Wp(rp) constraints in order to es- 
timate the HOD parameters and their errors by exploring the full 
parameter space of the model. 

In our HOD model the occupation number is parameterized 



(A^gall'w) = <A^cenl'«)(l + <A^satl'«)) 



(13) 



where (A^cenl'w) and (A^satl'w) are the average number of central 
and satellite galaxies in a halo of mass m. This model explicitly 
assumes that the first galaxy in haloes, when haloes have reached 
a sufficient mass, has to be central. Central and satellite galaxy 
occupations are defined as in (Zheng et al. 2005): 



(A^cenl'w) = ^ 



1 + erf 

m — Mo 
Ml 



log m - log Mn 



(14) 
(15) 



where M,nin, criog Mq, Mi , and a are the HOD parameters. The 
parameter Mo is generally poorly constrained and we decided in 
this analysis to fix Mo = Mj^in (see also White et al. 201 1; de la 
Torre & Guzzo 2012). 

In the HOD framework, the galaxy power spectrum or two- 
point correlation function can be written as the sum of two com- 
ponent: the 1-halo term that describes the correlations of galax- 
ies inside haloes and the 2-halo term that characterises the cor- 
relations of galaxies sitting in different halos (e.g. Abbas et al. 
2010). We follow the formalism of van den Bosch et al. (2012) 
to define the projected correlation in the context of this model. In 
particular we use their improved prescriptions for the treatment 
of the halo-exclusion and residual redshift-space distortions ef- 
fects on Wp(rp), induced by the finite n^n^^ values used in the data 
(van den Bosch et al. 2012). We use the halo bias factor and mass 
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Fig. 13. Evolution of the B-band absolute magnitude-dependent HOD. 
The curves show {Ng.^](m\z, Mg)) for values of x = Mg - 51og(/i) + z 
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function of Tinker et al. (2008) and Tinker et al. (2010) respec- 
tively, and assume that satellite galaxies trace the mass distribu- 
tion within haloes. We make the assumption of a NFW (NavaiTO 
et al. 1996) radial density profile and use the concentration-mass 
relation obtained by Prada et al. (2012) from the MultiDark sim- 
ulation. The details of the implementation of the HOD model 
are given in de la Torre at al. (in preparation). 

We present in Fig. 11 the measurements and best-fitting 
HOD models for the four different volume-limited absolute 
magnitude-threshold samples. We find that the model repro- 
duces the observations well. To have a global characterisation 
of the clustering properties of galaxies in VIPERS, we extend 
this modelling to two additional redshift bins at 0.5 < z < 0.7 
and 0.9 < z < 1.1. The best-fitting Mj^in and Mi parameters for 
the different sub-samples are shown in Fig. 12 and compared to 
previous measurement in the same range of redshift and number 
density. Because in the different surveys the subsamples are not 
selected with the same absolute magnitude band of selection, it 
is convenient to compare the HOD parameters in terms of red- 
shift and the number density probed by each sample. Note that 
here we compare measurements only from analyses using the 
same HOD parameterization, although the exact implementation 
of the models can differ slightly. The VIPERS sample allows us 
to constraint these parameters with an unprecedented accuracy 
over the redshift range 0.5 < z < 1.1. Our results are consis- 
tent with previous measurements, in particular with the DEEP2 
(Zheng et al. 2007) and CFHTLS (Coupon et al. 2012) analyses. 
Our HOD analysis is aimed at modelling the global clustering 
properties in VIPERS, but we refer the reader to Marulli et al. 
(2013) and de la Torre et al. (2013, in preparation) for detailed 
analysis and interpretation of the luminosity and stellar depen- 
dence of galaxy clustering and luminosity-dependent halo occu- 
pation respectively. 

We use the derived HOD parameters to define a global 
luminosity- and redshift-dependent occupation number which is 
then used to create accurate HOD mocks of the survey. In order 
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to interpolate between the different redshifts we assume a global 
luminosity evolution proportional to redshift, so that the magni- 
tude threshold values scale linearly with redshift (Brown et al. 
2008; Coupon et al. 2012). We find that one can approximate 
<A^gai(m|z, Mb)) using Eq. 13 with 



logMniinW = 10.61 exp(1.49 



-24.66-a: 



0-iogmW = 0.06exp(-0.08jc + 0.34) 
Mo(x) = M^Jx) 
Mi(x) = 13.5M„i„W 
a(x) = 0.29exp(-0.05jt; + 0.38) , 



(16) 
(17) 
(18) 
(19) 
(20) 



where x = Mg - 5 log{h) + z- M^in and Mi are found to be 
strongly correlated in such a way that Mj is approximately equal 
to 10-20 times M„im depending on the redshift probed and the 
model implementation (e.g. Beutler et al. 2013). In our analysis 
we find that Mi{x) can be approximated by 13.5 times Mj^mix). 
The function {Ngai(m\z, Mb)} is shown in Fig. 13 for the different 
values of x probed with VIPERS. We checked the consistency 
of this parameterization and verify that the Wp(rp) predicted by 
the mocks and the that measured are in good agreement for all 
probed redshift and luminosity thresholds. 

7. Redshift-space distortions 

The main goal of VIPERS is to provide with the final sample 
accurate measurements of the growth rate of structure in two 
redshift bins between z = 0.5 and z = 1.2. The growth rate of 
structure / can be measured from the anisotropics observed in 
redshift space in the galaxy correlation function or power spec- 
trum. Although this measurement is degenerate with galaxy bias, 
the combination /erg is measurable and still allows a fundamen- 
tal test of modifications of gravity since it is a mixture of the 
differential and integral growth. In this Section, we present an 
initial measurement of /erg from the VIPERS first data release. 

7.1. Method 

With the first epoch VIPERS data we can reliably probe scales 
below about 35/? ' Mpc. The use of the smallest non-linear 
scales, i.e. typically below 10 Mpc, is however difficult be- 
cause of the limitations of current redshift-space distortion mod- 
els, which cannot describe the non-linear effects that relate the 
evolution of density and velocity perturbations. However, with 
the recent developments in perturbation theory and non-linear 
models for RSD (e.g. Taruya et al. 2010; Reid & White 2011; 
Seljak & McDonald 2011), we can push our analysis well into 
mildly non-linear scales and obtain unbiased measurements of 
/cTg while considering minimum scales of 5 - 10 /z ' Mpc (de la 
Torre & Guzzo 2012). 

With the VIPERS first data release, we perform an initial 
redshift-space distortion analysis, considering a single redshift 
interval of 0.7 < z < 1.2. We select all galaxies above the mag- 
nitude limit of the survey in that interval. The effective pair- 
weighted mean redshift of the subsample is z = 0.8. The mea- 
sured anisotropic correlation function ^{rp,n) is shown in the 
top panel of Fig. 14. We have used here a linear binning of 
Arp - ts.K - Ihr^ Mpc. One can see in this figure the two main 
redshift-space distortion effects: the elongation along the line- 
of-sight, or Finger-of-God effect, which is due to galaxy ran- 
dom motions within virialized objects and the squashing effect 
on large scales, or Kaiser effect, which represents the coherent 
large-scale motions of galaxies towards overdensities. The latter 
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Fig. 14. Anisotropic correlation functions of galaxies at 0.7 < z < 
1.2. The top panel shows the results for the VIPERS first data release, 
deduced by the Landy-Szalay estimator counting pairs in cells of side 
1 /i"' Mpc. The lower two panels show the results of two simulations, 
which span the 68% confidence range on the fitted value of the large- 
scale flattening (see Section 7.4). 
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Fig. 15. Covariance matrix for the redshift-space distortions analy- 
sis. While the lower triangle represents the covariance matrix obtained 
directly from the MD mocks, the upper triangle corresponds to that ob- 
tained after performing shrinkage estimation. 

effect is the one we are interested in since its amplitude is di- 
rectly related to the growth rate of perturbations. Compared to 
the first measurement at such high redshift done with the VVDS 
survey (Guzzo et al. 2008), this signature is detected with im- 
pressive signal-to-noise, with the flattening being apparent to 
rp > 30/!"' Mpc. 

The two-dimensional anisotropic coiTelation has been exten- 
sively used in the literature to measure the growth-rate param- 
eter. However, with the increasing size and statistical power 
of redshift surveys, an alternative approach has grown in im- 
portance: the use of the multipole moments of the anisotropic 
coiTelation function. This approach has the main advantage of 
reducing the number of observables, compressing the cosmolog- 
ical information contained in the correlation function. In turn, 
this eases the estimation of the covariance matrices associated 
with the data. We adopt this methodology in this analysis and fit 
for the two first non-null moments ^o{s) and ^2(s), where most 
of the relevant information is contained, and ignore the contri- 
butions of the more noisy subsequent orders. The multipoles 
moments are measured from ^is,fi) which is obtained exactly 
as for ^(rp,7T), except that the redshift-space separation vector 
s is now decomposed into the polar coordinates (s,fi) such that 
rp = s{l - ji^yi^ and n - sfi. The multipole moments are related 
to ^is,fi) as, 

^eis) = ^— J Jis,^l)Lci^l)d^i, (2i) 

where L( is the Legendre polynomial of order i. In practice the 
integration of Eq. 21 is approximated by a Riemann sum over the 
binned ^{s,fj.). We use a logarithmic binning in s of Alog(i) = 
0.1 and linear binning in fi with Ajj - 0.02. 

7.2. Covariance matrix, error estimation, and fitting 
procedure 

The different bins in the observed correlation function and as- 
sociated multipole moments are correlated to some degree, and 
this must be allowed for in order to fit the measurements with 



theoretical models. We estimate the covariance matrix of the 
monopole and quadrupole signal using the MultiDark (MD) and 
Pinocchio (PN) HOD mocks. The genetic elements of the matrix 
can be evaluated as 

= Z {y'^''^ - y^'^^) {y'^'i^ - y^'^^) ^^2) 

where Nr is the number of mock realizations, y{s) is the quantity 
of interest, and the indices /, j run over the data points. 

The number of degrees of freedom in the multipole moments 
varies between 11 and 15 depending on the scales considered. 
Because we have only 26 MD mock realizations, the covariance 
matrix elements cannot be constrained accurately with the MD 
mocks only: the covariance matrix is unbiased, but it can have 
substantial noise. In order to mitigate the noise and obtain an ac- 
curate estimate of the covariance matrix, we apply the shrinkage 
method (Pope & Szapudi 2008), using the covariance matrix ob- 
tain with the 200 PN mocks as the target matiix. The PN mocks 
are more numerous and therefore each element of the associated 
covariance matrix is very well constrained, although the covari- 
ance may be biased to some extent. This bias is related to inaccu- 
racies in the predicted moments, which are mainly driven by the 
limited accuracy of the Zel'dovich approximation used in the PN 
mocks to predict the peculiar velocity field. The shrinkage tech- 
nique allows the optimal combination of an empirical estimate 
of the covariance with a target covariance, minimising the total 
mean squared error compared to the true underlying covariance. 
An optimal covariance matrix C is then obtained with 

C^AT + (l-A)S, (23) 

where A is the shrinkage intensity and the target T and empirical 
S covariance matrices correspond respectively to those obtained 
from the PN and MD mocks. A is calculated from (Pope & Sza- 
pudi 2008) 

^ ^ E,.;Cov(5,7,^,7)-Cov(r,v,^,v) ^^^^ 

where Cov(A,j, B,y) stands for the covariance between the ele- 
ments of the matrices A and B. We note that, because 
the empirical and target matrices are independent, the term 
Co\{Tij,S ij) vanishes in the numerator of Eq. 24. The effect 
of shrinkage estimation on the MD covaiiance matrix is shown 
in Fig. 15. 

To measure the growth rate of structure we perform a max- 
imum likelihood analysis of the data given models of redshift- 
space distortions by adopting the likelihood function X.: 

N„ N„ 

-21n^ = ZZ^'^'V^' ^25) 

where A^^ is the number of points in the fit, A is the data-model 
difference vector, and C is the covariance matrix. The likelihood 
is performed on the quantity y{s) - s^^ds), rather than simply 
y{s) - ^[{s) to reduce the range of variations of multipole values 
at different s in the fit. In the end, the quantity which is matched 
with model predictions is the concatenation of s^^q and for 
the set of separations considered. 

As a final remark, we note that we use the direct inverse of 
the covaiiance matrix without applying the correction discussed 
by Hartlap et al. (2007), as it is not clear how the size of the cor- 
rection is affected by the shrinkage estimation technique. The 
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Fig. 16. Systematic and statistical errors obtained from the mock 
samples, for different models and minimum scales s„,j„ in the fit. The 
horizontal dashed line shows the expectation value while the shaded 
area marks the 15% region around the fiducial value. Note that the 
points have been slightly displaced horizontally around j„„„ values of 
5 Mpc, 6 A"' Mpc, and 7 A"' Mpc to improve the clarity of the figure. 



resulting errors derived from the likelihood are well matched to 
the distribution of best-fit values from the mocks, which gives us 
confidence that only a small correction, if any, would be neces- 
sary. 

7.3. Models 

The formalism that describes redshift-space anisotropics in the 
power spectrum can be derived from writing the mass density 
conservation in real and redshift space (Kaiser 1987). In par- 
ticular, in the plane-parallel approximation, which is assumed 
in this analysis, the anisotropic power spectrum of mass has the 
general compact form (Scoccimarro et al. 1999) 

[6(x) + v^f0(x)mx') + v^fe(x')]) (26) 

where v - k\\lk, u\\{v) - -v\\(j-)l{faH{a)), V||(r) is the line-of- 
sight component of the peculiar velocity, 6 is the density field, 
9 is the divergence of the velocity field, Awy - u\\{x) - u\\{x') 
and r = X - x' . Although exact, Eq. 26 is impractical for direct 
use on redshift survey measurements and several models have 
been proposed to approximate it. In the assumption that galaxies 
linearly trace the underlying mass density field with a bias b, we 
can build three empirical models. These take the form. 



Pl(k,v)^D(kvcr,)PK(k,v;f,b) 
where, 

f exp{-{kvo-yf) 



(27) 



D{kvo-y) = 



and, 

PK(k,yJ,b) = 

b^Pss(k) + 2v^fbPss(k) + v^fPss(k) (model A) 



b'Pssik) + 2v\fbPse{k) + v^fPgeik) (model B) 



b'Psdk) + 2v\fbPsgik) + vy'Pgg(k) 

+CA(k, v;f,b) + CB(k,v;f,b). 



(model C) 



In these equations Pgg, Pgg, Pgg are respectively the non-linear 
mass density-density, density-velocity divergence, and velocity 
divergence-velocity divergence power spectra and cr,, is an ef- 
fective pairwise velocity dispersion that we can fit for and then 
treat as a nuisance parameter. The expressions for Ca(A;, v; /, b) 
and Csik, v; /, b) are given in appendix A of de la Torre & 
Guzzo (2012). These empirical models can be seen in config- 
uration space, as a convolution of a damping function Dikficr^,), 
which we assume to be Gaussian or Lorentzian in Fourier space, 
and a term involving the density and velocity divergence cor- 
relation functions and their spherical Bessel transforms. While 
the first term essentially (but not only) describes the Finger-of- 
God effect, the second, Pxik, v, b), describes the Kaiser effect. 
We note that model A is the classical dispersion model (Fisher 
et al. 1994; Peacock & Dodds 1994) based on the linear Kaiser 
(1987) model; model B is the generalisation proposed by Scocci- 
marro (2004) that accounts for the non-linear coupling between 
the density and velocity fields, making explicitly appearing the 
velocity divergence auto-power spectrum and density-velocity 
divergence cross-power spectrum; model C is an extension of 
model B that contains the two additional correction terms pro- 
posed by Taruya et al. (2010) to correctly account for the cou- 
pling between the Kaiser and damping terms. We refer the reader 
to de la Torre & Guzzo (2012) for a thorough description of these 
models. 

In the end, the model ^f{s) are obtained from their Fourier 
counterparts. 



2^ + 1 

P^(k) = ^- J ^dvPl(k,v)Ldy), 



dk 
2n^ 



eP',{k)jt{ks\ 



(28) 



(29) 



l/(l+(fcvcr,)2) 



where j[ denotes the spherical Bessel functions. 

The redshift-space distortion models involve the knowledge 
of the underlying mass non-linear power spectra of density and 
velocity divergence at the effective redshift of the sample. Al- 
though the real-space non-linear correlation function of galaxies 
can be recovered from the deprojection of the observed projected 
correlation function (e.g. Bianchi et al. 2012) and thus be used 
to some extent as an input of the model (e.g. Hamilton 1992), 
it is not feasible for the more advanced models which involve 
the velocity divergence power spectrum. The non-linear power 
spectra can however be predicted from perturbation theory or 
simulations for different cosmological models. In the case of a 
ACDM cosmology, the shape of the non-linear power spectra 
depends on the parameters P - {Q.„i,Q.h,h,ns,cri) and can be 
obtained to a great accuracy from semi-analytical prescriptions 
such as HALOFIT. In this analysis we use the latest calibration 
of HALOFIT by Bird et al. (2012) to obtain Pss and use the fit- 
ting functions of Jennings (2012) to predicted Pgg and Pgg from 
Pss- The latter fitting functions are accurate at the few percent 
level up to A: ^ 0.3 at z = 1 . 
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Fig. 17. Monopole and quadrupole moments of the redshift-space 
correlations, as a function of scale. The shallow curves show the results 
for the 26 individual MultiDark simulation mocks; the points are for the 
measured VIPERS data at 0.7 < z < 1.2, with assigned error bars based 
on the scatter in the mocks. The solid and dotted lines correspond to the 
best fitting models to the data for model B with Gaussian or Lorentzian 
damping function respectively. 

In the models, the bias and growth rate parameters b and 
/ are degenerate with the normalization of the power spectra 
parameter erg. Thus, in practice only the combination of bo-a 
and /cTg can be constrained if no assumption is done on the ac- 
tual value of (Tg. This can be done by renormalizing the power 
spectra in the models so that z) — » P.ix{k,z)/cr^(z), thus 

reparameterizing the models such that the parameters (b, f) are 
replaced by {bcr^,fcr^) in Eq. 27. This parameterization can 
be used for model A and B, although it is not for model C in 
which the correction term Ca involves the additional combi- 
nations: b^fcr'^, bf^cr^, and f^cr'^ (see Taruya et al. 2010; de 
la Torre & Guzzo 2012). The correction term Ca, which par- 
tially describes the effects of the non-linear coupling between 
the damping and Kaiser terms, mostly affects the monopole and 
quadrupole moments of the redshift-space power spectrum on 
scales of > 0.1 (Taruya et al. 2010). Therefore, in principle Ca 
could help breaking the degeneracy between / and cr^ although 
this has to be verified in detail. In the end, in the case of model 
C we decided to treat (/, b, as, cr,.) as separate parameters in the 
fit. 

7.4. Detailed tests against mock data 

We perform the redshift-space distortion analysis of the VIPERS 
data in the context of a flat ACDM cosmological model. Be- 
fore considering the redshift-space distortions in the data, we 
first test the methodology and expected errors on /erg using the 
mock samples. We fix the shape of the mass non-linear power 
spectrum to that of the simulation (since the observed real-space 
correlations are of high accuracy) and perform a likelihood anal- 
ysis of each individual MD mock. In the case of model C we also 
fix the normalisation of the power spectrum as discussed above. 
The distribution of best-fitting /erg values gives us a direct esti- 
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Fig. 18. Marginalized likelihood distribution of fcrg in the data (solid 
curve) and distribution of fitted values of fag for the 26 individual Mul- 
tiDark simulation mocks (histogram). These curves show a preferred 
value and a dispersion in the data that is consistent at the Icr level with 
the distribution over the mocks. 



mate of the probability distribution function of the parameter for 
a given fitting method, and serves as a check on the errors from 
the full likelihood function. We estimate the median and 68% 
confidence region of the distribution. These are shown in Fig. 
16 for the different models presented in the previous section and 
for various minimum scales 5„„„ in the fit. 

Model A is known to be the most biased model (e.g. Oku- 
mura & Jing 2011; Bianchi et al. 2012; de la Torre & Guzzo 
2012) and our results confirm these findings. We thus decided 
not to describe in the following the detailed behaviour of this 
model and focus on models B and C. We find that in general 
model B tends to be less biased than model C, which is surpris- 
ing at first sight as model C is the most advanced and supposed 
to be the most accurate (Kwan et al. 2012; de la Torre & Guzzo 
2012). This could be due to the quite restricted scales that we 
consider and the limited validity of its implementation on scales 
below s ^ 10 /i"' Mpc, as the maximum wavenumber to which 
we can predict Pgg and Pge is about k = 0.3. We defer the investi- 
gation of this issue to the redshift-space distortion analysis of the 
final sample and concentrate here on model B. The shape of the 
damping function in the models also affects the recovered /erg, 
as expected given the minimum scales we consider, although in 
the case of model B the change in /erg is at most 5%. Includ- 
ing smaller scales in the fit reduces the statistical error but at the 
price of slightly larger systematic error Therefore from this test 
we decided to use model and a compromise value for the mini- 
mum scale of s,„i„ - 6h^^ Mpc. 

7.5. Tine VIPERS result for the growth rate 

These comprehensive tests of our methodology give us con- 
fidence that we can now proceed to the analysis of the real 
VIPERS data and expect to achieve results for the growth rate 
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Fig. 19. A plot of /cTg versus redshift, showing VIPERS result contrasted with a compilation of recent measurements. The previous results from 
2dFGRS (Hawkins et al. 2003), 2SLAQ (Ross et al. 2007), VVDS (Guzzo et al. 2008), SDSS LRG (Cabre & Gaztafiaga 2009; Samushia et al. 
201 1), WiggleZ (Blake et al. 201 1), BOSS (Reid et al. 2012), and 6dFGS (Beutler et al. 2012) surveys are shown with the different symbols (see 
inset). The solid curve corresponds to the prediction for General Relativity in a ACDM model with WMAP9 parameters, while the dashed, dotted, 
and dot-dashed curves are respectively Dvali-Gabadaze-Porrati (Dvali et al. 2000), f{R), and coupled dark energy model expectations. For these 
models, the analytical growth rate predictions given in di Porto et al. (2012) have been used. 



that are robust, and which can be used as a trustworthy test of 
the nature of gravity at high redshifts. 

As explained earlier, we assume a fixed shape of the mass 
power spectrum consistent with the cosmological parameters ob- 
tained from WMAP9 (Hinshaw et al. 2012) and perform a max- 
imum likelihood analysis on the data, considering variations in 
the parameters that are not well determined externally. The best- 
fitting models are shown in Fig. 17 when considering either a 
Gaussian or a Lorentzian damping function. Although the mock 
samples tend to slightly prefer models with Lorentzian damping 
as seen in Fig. 16, we find that the Gaussian damping provides 
a much better fit to the real data and we decided to quote the 
corresponding /erg as our final measurement. 

We measure a value of 

fiz = 0.8)o-8(z = 0.8) = 0.47 + 0.08, (30) 

which is consistent with the General Relativity prediction in a 
flat ACDM Universe with cosmological parameters given by 
WMAP9, for which the expected value is /(0.8)o-8(0.8) = 0.45. 
The marginalised likelihood distribution of /erg is shown super- 
imposed on the mock results in Fig. 18. We see that the preferred 
values of the growth rate are consistent with the mocks, in terms 
of the width of the likelihood function being comparable to the 
scatter in mock fitted values. To illustrate the degree of flatten- 
ing of the anisotropic correlation function induced by structure 
growth, we show in the middle and bottom panels of Fig. 14 
^{rp, n) for two MD mocks for which the measured fcr^ roughly 
coincide with the Icr limits around the best-fit /erg value ob- 
tained in the data. We therefore conclude that the initial VIPERS 
data prefer a growth rate that is fully consistent with predictions 
based on standard gravity. Our measurement of /erg is also in 
good agreement with previous measurements at lower redshifts 



as shown in Fig. 19. In particular, it is compatible within Icr with 
the results obtained in the VVDS (Guzzo et al. 2008) and Wig- 
gleZ (Blake et al. 2011) surveys at a similar redshift, although 
WiggleZ measurements tend to suggest /erg values diminishing 
with redshift faster than expected in standard gravity (but see 
Contreras et al. 2013). 



8. Conclusions 

We have analysed in this paper the global real- and redshift- 
space clustering properties of galaxies in the VIPERS survey 
first data release. We have presented the selection function of 
the survey and the corrections that are needed in order to derive 
estimates of galaxy clustering that are free of observational bi- 
ases. This has been achieved by using a large set of simulated 
mock realizations of the survey to quantify in detail the system- 
atics and uncertainties on our clustering measurements. 

The first data release of about 54000 galaxies at 0.5 < z< 1.2 
in the VIPERS survey allows a measurement of the real-space 
clustering of galaxies through the measurement of the projected 
two-point coiTelation function, to an unprecedented accuracy 
over 0.5 < z < 1.2. This permits detailed modelling of the 
halo occupation distribution at these redshifts to be carried out. 
From an initial HOD modelling of B-band luminosity selected 
samples, we have been able to accurately determine the charac- 
teristic halo masses for halo occupation in the redshift interval 
0.5 < z < 1.0. These measurements are invaluable for creating 
realistic synthetic mock samples. 

The main goal of VIPERS is to provide an accurate mea- 
surement of the growth rate of structure through the character- 
isation of the redshift-space distortions in the galaxy clustering 
pattern. With the first data release we have been able to provide 
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an initial measurement of fcr% at z = 0.8. We find a value of 
/(Tg = 0.47 ± 0.08 which is in agreement with previous mea- 
surements at lower redshifts. This allows us to put a new con- 
straint on gravity at the epoch when the Universe was almost half 
its present age. Our measurement of /erg is statistically consis- 
tent with a Universe where the gravitational interactions between 
structures on 10 /z"^ Mpc scales can be described by Einstein's 
theory of gravity. 

The present dataset represents the half-way stage of the 
VIPERS project, and the final survey will be large enough to 
subdivide our measurements and follow the evolution of fa% out 
to redshift one. This will allow us to address some issues such 
as the suggestion from the WiggleZ measurements that /erg falls 
faster than expected at z > 0.5. Our measurement at z = 0.8 
already argues against such a trend to some extent, but the larger 
redshift baseUne and tighter errors from the final VIPERS dataset 
can be expected to deliver a definitive verdict on the the high- 
redshift evolution of the strength of gravity. 
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Appendix A: Galaxy mock catalogue construction 

We provide in this appendix some details about the method that 
we used to create realistic galaxy catalogues based on the Halo 
Occupation Distribution (HOD) and Stellar-to-Halo Mass Rela- 
tion (SHMR) formalisms. From the MultiDark simulation and 
Pinocchio halo lightcones described in Section 5.1, we created 



two types of galaxy mock catalogues: one containing B-band 
absolute magnitudes and associated quantities, and a second one 
containing stellar masses. We note that the stellar mass mock 
catalogues have not been explicitly used in this analysis, but in 
the accompanying VIPERS analyses of Marulli et al. (2013) and 
Davidzon et al. (2013). 

For the first set of catalogues we use the HOD formalism and 
populated dark matter haloes according to their mass by specify- 
ing the absolute B-band magnitude-dependent halo occupation. 
We parametrised the latter using Eq. 13 and used the HOD pa- 
rameters obtained from the data and given in Section 6. We po- 
sitioned central galaxies at halo centres with probability given 
by a Bernoulli distribution function with mean taken from Eq. 
14 and assigned host halo mean velocities to these galaxies. The 
number of satellite galaxies per halo is set to follow a Poisson 
distribution with mean given by Eq. 15. We assumed that satel- 
lite galaxies follow the spatial and velocity distribution of mass 
and randomly distributed their halo-centric radial position so as 
to reproduce a Navarro et al. (1996) (NFW) radial profile, 



PNFw(r\m) oc 



Cdm(m)r 
rvim) 



1 + 



Cd,n(m)l' 

rvim) 



(A.1) 



where c^m is the concentration parameter and ry{m) is the virial 
radius defined as 



ry{m) = 



3m 



47rpiz)ANL 



1/3 



(A.2) 



In this equation, p(z) is the mean matter density at redshift z and 
Anl = 200 is the critical overdensity for virialisation in our def- 
inition. We assumed the mass-concentration relation of Bullock 
et al. (2001): 



Cdm(m,Z) 



Co / _m y 
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(A.3) 



where co = 11 and m, is the non-linear mass scale at z = de- 
fined such as cr(m,, 0) = 6c. Here 6c and cr(m, 0) are respectively 
the critical overdensity (we fixed 6c = 1.686) and the standard 
deviation of mass fluctuations at z = 0. The latter is defined as 



cP-{m,z) = 

Jo 



dkk^P(k,z) 



(A.4) 



where R = [3m/ (4;rp(z))]'^^, P(k,z) is the linear mass power 
spectrum at redshift z in the adopted cosmology, and Wix) is the 
Fourier transform of a top-hat filter. In order to assign satellite 
galaxy velocities, we assumed halo isotropy and sphericity, and 
drew velocities from Gaussian distribution fimctions along each 
Cartesian dimension with velocity dispersion given by (van den 
Bosch et al. 2004): 
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(A.5) 
(A.6) 



where t//(r) is the gravitational potential, G is the gravitational 
constant, fix) = ln(l + x) - x/il + x), and 



I(x)= f 

Jx 



f{t)dt 



(A.7) 



In these mocks, the absolute Z?-band magnitude for each 
galaxy was obtained following Skibba et al. (2006). From the 
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mean rest-frame B — i' colour and K-corrections observed in the 
data we then derived absolute and apparent j'-band magnitudes 
for each simulated galaxy. 

For the mock catalogues with stellar masses, we followed 
the SHMR approach which is based on the assumption of a 
monotonic relation between halo/subhalo masses and the stellar 
masses of the galaxies associated with them. We first popidated 
the haloes in the lightcones with subhaloes. For this we ran- 
domly distributed subhaloes around each distinct halo following 
a NFW profile so that their number density satisfies the subhalo 
mass function of GiocoU et al. (2010): 



rflnmsub 



= ^orexp(-j8r), (A.8) 



where m^a\, is the subhalo mass, ^ = mf^ut^/m, a = -0.8, /3 - 
12.2715, and A^o = 0.18. We then assigned a galaxy to each halo 
and subhalo, with a stellar mass given by the SHMR of Moster 
et al. (2013). The galaxy velocities were assigned in a similar 
way as for the HOD catalogues, with galaxies associated with 
distinct haloes and subhaloes being considered as central and 
satellite galaxies respectively. 
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